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In these notes we present preliminary results on quantum-like algorithms where tensor product is 
replaced by geometric product. Such algorithms possess the essential properties typical of quantum 
computation (entanglement, parallelism) but employ additional algebraic structures typical of geo- 
metric algebra — structures absent in standard quantum computation. As a test we reformulate in 
Geometric Algebra terms the Deutsch-Jozsa problem. 



Quantum algorithms are based on tensor products. Common wisdom states that tensor products are needed for 
entanglement. A similar situation was encountered in the early 1990s in connectionist systems, and led to the paradigm 
known as tensor product representations However, nowadays the cognitive science community seems to depert 
from tensor product representations in favor of their "compressed forms" such as Binary Spatter Codes (BSC) |2j or 
Holographic Reduced Representations (HRRs) 0. 

The main reason why BSC and HRRs replace tensor product representations is that tensor multiplication expands 
dimensions of the associated liner spaces (tensor product of two n-tuples is an n 2 -tuple). HRRs, for example, replace 
tensor product with circular convolution, an operation that does not change the dimension (circular convolution of two 
n-tuples is again an ri-tuple). Circular convolution is often referred to as a compressed form of the tensor product. 
Similarly, in BSC one replaces tensor products by XORs of binary strings. Quite recently, following the general 
program of investigating similarities and differences between AL semantic analysis, and quantum information jj], we 
have reformulated BSC in terms of Geometric Algebra (GA) [HuIl]]- This reformulation was made possible by the 
observation that XOR has a natural representation at the level of geometric product. 

In these notes we present a similar reformulation of the Deutsch-Jozsa algorithm 8] . As one can see, after minor 
modifications the GA algorithm works analogously to the quantum one. Accordingly, it is possible that GA algorithms 
can perform more general tasks until now reserved for quantum computation. The fact that it was easy to reformulate 
in a GA way the Deutsch-Jozsa problem is very encourageing. 

Finally, let us mention that certain attempts of using GA for the purposes of quantum computation can be found 
in the literature seems that the approaches discussed so far reduce GA to the level of unitary 

operations or density matrices, i.e. objects that have a natural operator representation. In our approach, even "pure 
states" are represented by elements of GA. This is why we can perform operations on pure states that have no 
couterpart in standard quantum computation. In this sense our geometric algorithm may be regarded as something 
conceptually in-between quantum algorithms and HRRs or BSC. 



I. INTRODUCTION 



II. ORIGINAL DEUTSCH-JOZSA ALGORITHM 



We assume there exists an oracle performing 



U f \x)\y) 



\x)\y®f(x)) 



(i) 



where f(x) £ {0, 1}. Now 
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If f(x)=0 then 
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If f(x) = l then 
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The two cases imply 

tf/|*)(|0)-|l>) = (-l)^k)(|0)-|l)) 

The Hadamard gate acts as follows 

U H \0) = ^(|0> + |1>) 
U H \1) = -^=(|0> - |1> 

Let 

U n+1 = U H ®---®U H 

v ' 

n+1 

Then 



C/ n+1 |0^_01) = -L= ]T ^•••^)(|o)-|i)) 

= 5>>(i°>-ii>) 

U f U n+1 \0_^01) = -jL= £ (-1)^1^(10) -|1)) 

1 2 "-l 

U n+1 U f U n+1 \ 0_^0l) = -p^t-l/W^Wjl) 

V ^ Ai...A n =0 

= ^ E E (-i) ,(Ai - A " ) (-i) Et ^» B «'iBi--B»>ii> 

Ai...A n =0 Si...S n =0 

= ^ E ( E (-i) /(Ai - a » ) (-i) e - a » b »)ibi...b„)|i> 

Bi...S n =0 A 1 ...A n =0 
(Bi...S n )^(Oi...O n ) Ai...A n =0 

+ E (-l^-^lOi.-.O^Il) 

A 1 ...A n =0 

= ••• + ^E(- 1 ) /(x) i°i--- «)i 1 ) 

£C = 

If /(x) = /(0) for any a; then 

RHS = ... + i r (-l)/(0) Y |0i...0 n )|l) = (-l)/(o)| 0l . ..0„)|1) 

x=0 
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If / is balanced then 

RHS = ••• + 4 £(-i) /(:c) |Oi...o n >|i> 



On /.^ 



E (-1: 



/(Ai...A„) 



(-l)E5UW)| fll ....B n >|l 



(S 1 ...B„)#(Oi...O„) Ai...A„=0 

It is sufficient to look at the |0 . . . 01) component to see if / is constant or balanced. 



(21) 
(22) 



III. GEOMETRIC ALGEBRA AND ITS BINARY PARAMETRIZATION 



Here and in the next section we repeat the presentation from [5|. 

Euclidean-space GA is constructed as follows. One takes an rt-dimcnsional linear space with orthonormal basis 
{ei, . . . , e n }. Directed subspaces are then associated with the set 



{Mi 



1 e n , ei2, ei3 . . . , e n -x,n-> •••> e i2...n}- 



(23) 



Here 1 corresponds to scalars, i.e. a 0-dimensional space. Then we have vectors (oriented segments), bivectors 
(oriented parallelograms), and so on. There exists a natural parametrization: 1 = eo...o, e± = eio...o, &2 = eoio...o, ■ ■ ■ , 
e i25 = e nooio...Oi . . . , ei2...n-i,n = e n...i> which shows that there is a one-to-one relation between an n-bit number 
and an element of GA. An element with k Is and n — k 0s is called a fc-blade. 

A geometric product of k 1-blades is a fc-blade. For example, ei248 = e\e2e±e%. Moreover, e n e m — —e m e n , if to 7^ n, 
and e n e n = 1, for any n. GA is a Clifford algebra 01 enriched by certain geometric interpretations and operations. 

Particularly interesting is the form of the geometric product that occurs in the binary parametrization. Let us work 
out a few examples: 



eiei 
ei2ei 

ei257e26 



eio...oeio...o = 1 = e o...o = e (io...o)e(io...o) 

eio...o e no...o = eieie 2 = e 2 = e io...o = e (io...o)e(no...o) 

eno...o e io...o — — ~^2&\ e \ = = — e io...o = — e (no...o)e(io...o) 



(24) 
(25) 
(26) 



enooioio...oeoioooio 



...0 = eie2e5e 7 e 2 e 6 = (-if l eie 2 e 2 e 5 e 7 e e = (-1) 2 '(-I) 1 eie2e 2 e 5 e 6 e 7 



(— I) 3 eie5e6e7 — (— l) 3 eiooomo...o — (^l) I?e (iiooioio...o)e(oioooio...o)- (27) 



The number D is the number of times a 1 from the right string had to "jump" over a 1 from the left one during the 
process of shifting the right string to the left. Symbolically the operation can be represented as 



01000100... 



11001010... 



(-1) J 



01000100. 
11001010. 



(-1) J 



The above observations, generalized to arbitrary strings of bits, yield 

1 BkAi 



1 \ ' n .../;. = (-1)£*<i fl » A «e 



Indeed, for two arbitrary strings of bits we have 

< — BiB 2 ■ ■ 



AlA 2 . . . A n 



(-1) 



(A, 



D 



01000100. . 

© 

11001010. . 

„)e(Bi...s n )< 



= (-1) D [ 10001110... 0] 



(28) 



3^2 

A X A 2 



■ B n 

■ A n 



where 



D = B X (A2 + ■ ■ ■ + An) + B 2 {A 3 + ■ ■ ■ + A n ) + ■ ■ ■ + B n -\A n = B k Ai 



(29) 



(30) 



k<l 
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IV. CARTAN REPRESENTATION 

In this section we give an explicit matrix representation of GA. We begin with Pauli's matrices 

<ri = ( 5 1 ) , ^ = (! _ o ! )' ct3 = ( J -°i ) • (31) 

GA of a plane is represented as follows: 1 = 2x2 unit matrix, e\ = a%, e-x = e 12 = <J\G<i = i(7 3 . Alternatively, we 
can write eoo = 1, em = <7i, eni = cr 2 , en = 1173, and 

/ aoo + iotix aio — iceoi 1 
a 00 eoo + ai ei + a ie i +anen = . (32) 

\ «10 + ia 01 a oo — ta ll 

This is equivalent to encoding 2 2 = 4 real numbers into two complex numbers. 

In 3-dimensional space we have 1 = 2x2 unit matrix, e\ = a±, e 2 = <7 2 , £3 = 03, &12 = = ^3, eis = (J1O3 = 
— i(72, e 2 3 = 0203 = «(7i, ei 2 3 = <7i<7 2 cr3 = i. 

Now the representation of 

E a _( «ooo + *am+ a 00 i + iano, at wo + ia n - ia ow - a W i \ , , 

\ CKinn + ICenn + iamn + aim . Cennn + idfin — anm — iai in / 



ABC=0.1 



is equivalent to encoding 2 3 = 8 real numbers into 4 complex numbers. 

An arbitrary n-bit record can be encoded into the matrix algebra known as Cartan's representation of Clifford 
argebras |T^ : 

e 2fc = (7i g> ■ ■ ■ ® gi g 2 (8) 1 (8) ■ ■ ■ ® 1 , (34) 

n — fe fc — 1 

ea/fe-i = pi • ■ ■ (J\ 03 1 — ■ 1 ■ (35) 

n — k k—1 

V. GA FORMULATION OF THE DEUTSCH-JOZSA PROBLEM 

Consider an (n+ l)-dimensional Euclidean space with orthonormal basis {ei, . . . e n _|_i}, and its associated GA. The 
basis vector e„+i in binary parametrization corresponds to en... 01. Recall that 

eA 1 ...A n+1 e Bl ...B n+1 = (-1) Ei <j B * Aj 'e (j 4 1 ... J 4„ +1 )e(Bi...Bn+i) ( 36 ) 

and, in particular, 



eAi...A„+i e o...oi — e J 4 1 ...A„,A„+i®i (37) 
eAi...A„ +1 eo...oio = (-l) A " +1 e J 4i...A„_i,A„©i,A„ + i (38) 



Consider 

1 

E n +i = E eA x ...A n+ i (39) 

Ai...A n +l=0 

1 

£'n+ieo...oio = E eAi..A„+i e o...oio (40) 
Ai...A„+i=0 
1 

E (-l)^ 1 ^...^.!,^©!,^ (41) 
A 1 ... J 4„ + 1 =0 

1 1 

E (-^ An+1 e Al ...A n ^,om,An +1 + E (-l) An+1 e Al ...A n - 1 ,i@i,An+i ( 42 ) 

Ai...A„_iA„ + i=0 Ai...A„_iA„ +1 =0 
1 1 

E (-l)^ n+1 e Al ...vi B _ 1 ,i,A n+1 + E (-l) An+1 eA 1 ...A„_ 1 ,o,A„ +1 (43) 

A 1 ...A n -iA n+1 =0 Ai...A„_ lJ 4„ + i=0 
1 

= E (-l) A " +1 ^...A„ +1 (44) 
Ai...A n+i =0 

The influence of £Vi+i on eo...oio is similar to (|11|) : 

l 

£„+ie ...oio = E ("l) A " +1 eA 1 ...A„ +1 (45) 

Ai...A„+i=0 
1 

= E (eA 1 ...A n o - e Al ...A n i) (46) 

Ai..A„=0 

i 1 

^lO^Ol) = -_ E (-l) A " +1 l^i---^+i> (47) 

V/! Ai...A n =0 
Now assume there exists an oracle Ef that performs 

E f e Al ...A n A n+1 = eA 1 ...A„,A n+ i®/(Ai...A„) = eAi...A„,A n+1 e ...o,/(A 1 ...A„) (49) 

Then 

i 



£/£ n+ ie ...oio = E E f( e A 1 ...A n o-eA 1 ...A n i) (50) 

Ai...A n =0 
1 

= E S / ( e Ai...A„0 - eAi...A„l) (51) 
Ai...A„=0 
1 

= E ( e Ai...A„,/(Ai...A„) - eAi...A„,-./(Ai...A„)J (52) 
Ai...A„=0 
1 

= E (-l) /(Al - A " ) (eA 1 ...A„o-eA 1 ...A„i) (53) 
Ai...A„=0 

In GA there exists an operation of reverse which reverses the order as follows: If X = e\e2 ■ ■ ■ ek-iek then the reverse 
of X is 

Xt = efc e fc _! . . . e 2ei = (-l) fc (^D/ 2 X (54) 
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By linearity we extend it to all multivectors. In binary parametrization the number k describes the number of Is in 
e Ai...A n+1 , i-e. k = Y^=\ Ar So consider 

1 

F n+1 = E e Ai...A„o ( 55 ) 

A-L...A n =0 

= E (-l) Kk - 1)/2 e Al ...A n0 (56) 

Ai...A n =0 

Here k — Y^j=i A? smce the last bit is 0. Now 

l 

F n+1 E f E n+1 e ... 010 = ]T (-l) f(Al - A - ] F n+1 (e Al ... An0 -e Al ... A ^ (57) 

A!...A„=0 
1 

- E (-l) /(Al - A " ) ^+ieA 1 ...^o + ... (58) 

Ai...A„=0 

= E (-l)'^"^ E (-l)* ( *- 1)/a e Bl ...B.oeA,..^o + ... (59) 

Ai...A„=0 Bi...B„=0 

= E (-l)^ 1 -^^ E (-l) fe(fe - 1)/2 (-l) 5: -' Afci3 ^(B,..B„o)e(^,..A„o) + ...(60) 

Ai...A„=0 Bi...B„=0 
1 

- E (-l) /(Al - A " ) eo...o + ... (61) 

Ax...A n =0 

The dots denote all those term where the binary indices contain at least one 1. The two powers of —1 have cancelled 
out since e Ai An+i e Al ...A n+1 = 1 = e ...o- Finally 

l 

F n+1 E f E n+1 e ... 010 = E (-l) /(Al - A " ) eo...o + --. (62) 

A!...A„=0 

Now let II project on 1 = e ...o- It follows that 

l l 
TrU(F n+1 E f E n+1 e .. mo ) = E (-l) f(Al - An) Tr e ...o = N E (-l)'^ 1 "^ (63) 

Ai...A n =0 Ai...A n =0 

( (-l)/(°-°)jV2 n if /is constant , . 

\ if / is balanced 

Here N = Tr 1 is the dimension of the representation of GA. We have achieved the same goal as the quantum 
algorithm. 

We have to point out at this moment a possible error one can make. Let us note that in the step 

Efe Al ,,, AnAn+1 = e Al ...A n ,A n+1 eo...oj(A 1 ...A n ) (65) 

we have Ef on the left and eo...oj{ Al ...A n ) on the right. It might appear that it would be simpler and more natural 
to write Ef on the right as well. However, this would be misleading since 

F n+1 Efe Al ...A n A n+1 = F n+ i(e Al ...A n A n+1 e ...o,f(A 1 ...A n )) (66) 
7^ (^F n+ ie Al ... AnAn+1 ^je j( Al An ) = EfF n+1 e Al ... AnAn+1 (67) 
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VI. EXPLICIT EXAMPLES 
A. Two bits 

GA of a plane consists of: 1 = 2x2 unit matrix, e\ = a\, e 2 = cr 2 , £\i = &ia 2 = Alternatively, we can write 



eoo = 1, eio = e i = cr 2 , en = 

( a o + iau aio - ioioi \ , ao ^ 

aooeoo + aioeio + «oieoi + anen = . . . (68) 

y a w + ictm «oo — ia n J 

*-(i:J!:{)- (69) 

F 2 = ej + el = e 00 + eio = 1 + o-i = ^ j j ^ • (70) 

E 2 e w = (e oa + eio + eoi + enj e w = e w + e 00 - en - e i (71) 

= (l + (Ji+ a 2 + io-s) cri = a\ + 1 - ia 3 - a 2 (72) 

E f E 2 e w = E f (e w + e 00 - e n - e i) (73) 

= e 1)0 ©/(i) + e ,oe/(o) - e iae/(i) - e o,ie/(o) (74) 

= e li/( i) + e 0>/ ( ) - e W (i) - e ^ /(0) (75) 



J. Case /(0) = /(l) = 



EfE 2 e w = e\j(i) + e ,/(o) - e i,-./(i) _ e o,-/(o) (76) 
= eio + e 00 - en - e i (77) 
= o-i + 1 - ia 3 - a 2 (78) 



1 - % 1 + i 
1-i 1 + i 



(79) 



TrlU^-E/^eio = 4 (81) 



2. Case /(O) = /(l) = 1 



E f E 2 e 10 = ei i/( i) + e 0j/(0) - ei^ /( i) - e ,-,/( ) (82) 

= en + eoi - eio - e 00 (83) 

Since this is minus the result from the previus subsection, we immediately get 

Tr YlF 2 E f E 2 e 10 = -4 (84) 
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3. Case /(O) = 0, /(l) = 1 



EfE 2 e w = e x + e 0i /( ) - e i,-./(i) - e o,-./(o) (85) 
= en + e o - e w - e i (86) 
= ia 3 + 1 - CTi - C7 2 (87) 



1+z -1 +i 
1-i 



(88) 



Alternatively 



fi^-iii . 1 1 + J j - 1 1 .V - SS) < 89 » 



Trn^/^eio = (90) 



F 2 E f E 2 e w = (l + CTi)(io- 3 + 1-0-1-0-2) (91) 

= (Z03 + 1 - Oi - 02) + (ioi03 + 01 - 1 - O1O2) (92) 

= (zo 3 + 1 - 01 - o 2 ) + (i(-«o 2 ) + 01 - 1 - ia 3 ) = (93) 

4. Case /(O) = 1, /(l) = 



EfE 2 e w = ei i/( i) + e ,/( ) - e w( i) - e ^ /(0) (94) 
= eio + e i - en - e 00 (95) 

This is minus the result from the previou section and therefore Tr UF 2 EfE 2 e w = 0. 

Summing up, constant functions were producing 2(— l)^ ^ 1 , and balanced functions implied 0, as it should be on 
general grounds. 

B. Three bits 

In 3-dimcnsional space we have 1 = 2x2 unit matrix, ei = 01, e 2 = 02, e3 = 03, ei2 = 01O2 = 203, ei 3 = 0103 = 
— io2, e23 = 0203 = iai, e\ 23 — o\a 2 a 3 = i. The operation Trll corresponds in this representation to taking the real 
part of trace (only e oo = 1 and e m = i have nonzero trace). 

Now the representation of a general element reads 

E_( c*ooo + + a ooi + iociw, aioo + iatou - iaoio - "101 \ /qg\ 
abc abc \ a wo + ia 01 i + ia ow + a W i, a oo + ioc\\i — aooi - iano / 



ABC=0,1 



£3- ^ ( 1 + z + z + l, 1 + z-l-z ) =2 ( 1 + z, J (97) 



ABC=0,1 



F 3 = e AB0 = e 000 + e 100 + e 010 + e 110 = 1 + °"l + ^2 + (oi0 2 ) f (98) 



AB=0,1 



1 — i 1 — i 

1 + 0-1 + o 2 + o 2 oi = 1 + 01 + o 2 - ia 3 = ( 1 ^ j ^ ) (99) 
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*--»(i:!S)(?^)-»(S!:0 « 

= ^e oo + eioo + e io + eooi + eno + e n + eioi +einje io (101) 

= eoooeoio + eiooeoio + eoio e oio + eooieoio + enoeoio + eoneoio + eioieoio + emeoio (102) 

= e i + eno + e 00 o - e n + eioo - e 00 i - e m - eioi (103) 

= cr 2 + a x a 2 + 1 - cr 2 er 3 + (Ti - <7 3 - o-i^fa - ^l*^ (104) 

= <t 2 + icr 3 + 1 - iffi + <7i - cr 3 - i + ia 2 (105) 
i + 1 - 1 - i, -i - i + 1 + 1 



+ -i + l + l-i 



(106) 



EfE 3 eoio — eoi^e/foi) + e ii,oe/(n) + e oo,oe/(oo) — e oi,ie/(oi) + eio,oe/(io) — e o,ie/(oo) — e ii4©/(n) ~ e io4®/(io) 
= e oi ,/(oi) + en,/(ii) + e oo,/(oo) — eoi,-,/(oi) + ei ,/(io) — e 00 ,^/(oo) — e n,^/(ii) — e io,^/(io) 
= e oo ,/(oo) - e o,-,/(oo) + e i ,/(oi) - e i,-,/(oi) + e io,/(io) - ei ,-,/(io) + e n,/(ii) - e n,-/(ii) (107) 



1. Case of constant f, /(00) = 



EfE 3 emo — eooo — eooi + eoio — eon + eioo — e ioi + eno — em (108) 

= 1 — <J3 + <J2 — io'i + <J\ + *o"2 + i&3 — i (109) 

-(1:^:1:1^)^(11) ("») 

F a E,E 3nm = (};;};;) 2(1-0 (J ;) an) 

-»c-o(l;J ;;:)(-) < ii2 » 

SRTrFaE/^seoio - 2(1 - i)2(l + *) = 8 = 2(-l)°2 2 (114) 
2. Case of constant f, /(00) = 1 

EfE 3 e aw — eooi — eooo + eon — eoio + eioi — eioo + em — eno (H5) 

MTr F 3 E f E 3 e 010 = -2(1 - i)2(l + t) = -8 = 2(-l) : 2 2 (116) 



5. Case o/ balanced f, /(00) = 0, /(10) = 



EfE 3 e aw — e o,/(oo) + eio,/(io) + e oi ,/(oi) + e n,/(n) — e oo^/(oo) — ei ,-,/(io) — e i,-,/(oi) — eu,-,f(u) 
= e oo + eioo + eon + em — e oi — eioi — e io — eno 
= e oo + eioo — e io — e oi + eon — eioi — eno + em 

::::^:;:*:i::Hi+»U:) <*»> 
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F 3 E f Ese 010 = 2(l + i)(° 2(1 + (118) 

SJTrJVE/.Eseoio = 3?2(1 + i)2(l + i) = K8i = (119) 
VII. FINAL REMARKS 

The above examples show that GA allows for a host of new mathematical tricks with respect to standard quantum 
computation. The representations of binary numbers are different. There is no distinction between "state vectors" 
and "operators" . One can multiply "state vectors" without increasing the dimension. In the above examples both 
2-bit and 3-bit problems were represented by 2 x 2 matrices, a fact showing that one may expect GA to involve less 
redundancy than standard tensor representations. One can speak of entanglement in GA representations even though 
the "states" are not tensored with one another. Here again one finds close analogies to what is known from HRRs 
and BSC. And, last but not least, it seems there is no general difficulty with translating quantum operations into GA 
forms, and one can expect all quantum algorithms to have GA analogues. 
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